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Spin transition in Gd 3 N@C 80 , detected by low-temperature on-chip SQUID 
technique 
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We present a magnetic study of the GdsN^Cso molecule, consisting of a Gd-trimer via a Nitrogen atom, 
encapsulated in a Cso cage. This molecular system can be an efficient contrast agent for Magnetic Resonance 
Imaging (MRI) applications. We used a low-temperature technique able to detect small magnetic signals 
by placing the sample in the vicinity of an on-chip SQUID. The technique implemented at NHMFL has 
the particularity to operate in high magnetic fields of up to 7 T. The GdaN^Cso shows a paramagnetic 
behavior and we find a spin transition of the GdaN structure at 1.2 K. We perform quantum mechanical 
simulations, which indicate that one of the Gd ions changes from a 8 5V/2 state (L = 0, S = 7/2) to a 7 F 6 
state (L = S = 3, J = 6), likely due to a charge transfer between the Cgo cage and the ion. 
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The metal-containing endohedral fullerenes molecules 
have drawn attention as they exhibit unusual material 
properties associated with charge transfer between the 
metal cores and the carbon cage. 1 The encapsulated 
metal cores are protected from the solvent or other agents 
by the carbon cage, which leads to potential biologi- 
cal applications. Specifically, the trigadolinium fullerene 
Gd3N@C8o (in Fig 1) has attracted wide interest, because 
it is identified as a safe candidate as contrast-enhancing 
agent in Magnetic Resonance Imaging; also, it can be im- 
planted with other ligands for medical applications 2 . The 
carbon cage separates the interior metal atoms from the 
environment, preventing their leakage and accumulation 
into the human body. Also, vibrational modes 3 indicate 
a bond formation between the cage and the GdaN. Stud- 
ies showed 4,5 that the GdaN cluster is bound to the Cso 
cage by a large binding energy of 13.63 eV which ensures 
that Gd atoms are not likely to break out of the cage. 
The three Gd atoms are bonded through their delocal- 
ized s and d states, with the s and p states of C atoms. 
This leaves the spin magnetic moment of the / electrons 
almost intact, with no appreciable spin density on the 
cage. This compound can improve the MRI relaxation 
enhancement by a factor of 20, leading to a new class of 
MRI tumor delineation techniques. 5 Such developments 
towards novel medical applications, require a fundamen- 
tal understanding of the electronic and magnetic proper- 
ties of these endohedral metalofullerenes. In this paper, 
we present a study of the spin properties of the GdaN 
core cluster by on-chip Superconducting Quantum Inter- 
ference Device (SQUID) measurement technique. 

Magnetic flux detectors using SQUID devices feature 
very high sensitivity and operate at low temperatures. 6 




FIG. 1. The sketch of the on-chip SQUID setup. A feedback 
coil is placed above an on-chip SQUID and the magnetic field 
is oriented parallel to it. The Gd3N@Cso sample (not at scale) 
is placed near the on-chip SQUID loop. The structure of the 
GdsN^Cso molecule is presented in the insert. 



Unlike commercial SQUIDs, where the sample is coupled 
to the detector via a pickup coil, the on-chip SQUID 
method couples the sample directly, by placing it on 
the same chip as the SQUID, for increased sensitivity 7 . 
The loop of a SQUID is interrupted by two Josephson 
junctions (see Fig. 1), which can be of planar-tunneling 
type 6 or nano-bridges 8 which are constrictions of a thin 
niobium film and are much thinner than tunneling junc- 
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tions. Thin films have a significantly larger in-plane criti- 
cal field, compared to the bulk value 9 . Consequently, the 
nano-bridge junctions are less affected by a large mag- 
netic field applied along the film plane. By aligning the 
magnetic field with high precision in the plane of a Nb 
SQUID, the on-chip SQUID technique employed here ex- 
tends the observation of quantum magnetic phenomena 
into a higher magnetic field range. Therefore, this im- 
proved technique is very promising in gathering more in- 
formation about effects such as spin state transitions and 
entanglements, the anisotropy-induced tunneling gap 10 , 
phonon bottleneck effects 11,12 and others. 

SQUID operation require a precise alignment of the 
magnetic field B. We can rotate B in SQUID's plane 
by combining fields generated by a 3- axis magnet B = 
B x + By + B z . Initially, the SQUID plane is roughly 
aligned along B z . We can tune the direction of the field 
along the crossing between the SQUID plane and the 
xz plane by applying a small amount of B x , which has 
been set roughly normal to the SQUID plane (similar 
approach applies for the yz plane, if using the B y com- 
ponent). The on-chip SQUID is made of a Nb film of only 
5.5 nm thickness ensuring a superconducting state even 
for 7 T of in-plane magnetic field 13 . The SQUID loop has 
a 2x2 fim 2 area and contains two parallel nano-bridge 
Josephson junctions. This technique allows us to place 
a small amount of GdaN^Cso particles in the vicinity of 
the SQUID loop for an increased magnetic coupling be- 
tween sample and SQUID. The entire device is thermally 
anchored to the mixing chamber of a dilution refrigerator 
to ensure an efficient cooling. 

The detection is done by analyzing the output volt- 
age following a series of current pulses, sent to the 
SQUID. If the current pulse is sufficiently strong, the 
SQUID switches into the normal state, generating de- 
tectable voltage pulses. Statistically, the ratio between 
the number of output voltage pulses and that of the ap- 
plied current pulses, gives P SW1 the switching probabil- 
ity 14 . SQUID switching into the normal state generates 
heat, which we minimize by operating the SQUID in 
pulse mode. In our experiment, current pulses are gener- 
ated by a pulsed voltage source with an in-series resistor. 
The switching probability P sw of the SQUID increases 
with the current amplitude flowing through the junc- 
tions. This can be observed in the contour plot of Fig. 2, 
showing sharp transitions between low and high switch- 
ing probabilities. We define the switching current I sw as 
the current corresponding to P sw = 50% and Fig. 2 shows 
a typical periodic modulation as a function of the mag- 
netic flux in the SQUID loop, with a period of $o = h/2e 
(h is Planck's constant and e the electron charge). 

Real time ( "on-the-fly" ) magnetic measurements are 
performed using a feedback mechanism, described below. 
An "working point" for the SQUID is chosen, for instance 
at the middle point in the modulation curve (see Fig. 2). 
During the measurement, a current pulse with a height 
corresponding to the working point is sent to the SQUID. 
It is important to note that the variation of the proba- 
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Feedback coil current (mA) 

FIG. 2. Contour plot showing P sw as a function of the feed- 
back current and the height of the switching pulse applied to 
SQUID. The sharp transition from low to high switching prob- 
ability are a characteristic of SQUID devices. The detector 
is kept at a working point (here P sw — 50%) by the feedback 
mechanism, relaying information on the compensated exter- 
nal flux, coming from a nearby sample. 



bility near the working point is very sharp, which is the 
key point behind the high flux sensitivity of a SQUID. 
We use a small superconducting coil placed above the 
SQUID (see Fig 1), to generate a flux compensating any 
flux changes caused by the magnetic sample. The feed- 
back coil current is proportionally related to the change 
in the magnetic moment of the sample. This method 
locks the SQUID at its working point and performs a 
feedback measurement. 

The described setup is used to measure the magnetiza- 
tion of Gd3N@C8o molecules, conglomerated in a small 
particle, while sweeping the applied field from negative to 
positive saturation ( — 1 to +1 T). The magnetic field was 
set at a sufficiently low sweeping rate to ensure a thermal 
equilibrium magnetization process. The high symmetry 
of the Cso cage does not favor the existence of magnetic 
anisotropy in the GdaN system, as confirmed by our re- 
versible magnetization curves (see below). The collective 
spin momentum of GdsN core cluster exhibit a paramag- 
netic behavior, as expected also from theoretical calcu- 
lations. 5 In the insert of Fig. 3 magnetization curves at 
different temperatures are presented. No hysteresis is ob- 
served for all temperatures, confirming that the sample is 
completely paramagnetic. The magnetization curves are 
flattened with the increase of temperature, which gives a 
smaller initial susceptibility x- To obtain the initial sus- 
ceptibility, we normalize the magnetization curves to the 
saturation value (at very low temperature and 1 T) and 
then take the derivative at zero field. A plot 1/x vs. T 
is shown in Fig. 3. We observe a linear dependence but 
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Temperature (K) 

FIG. 3. The inverse of the initial susceptibility (black dots) as 
a function of temperature, calculated from normalized magne- 
tization curves (see insert). The continuous line is a fit based 
on Eqs. 1 and 2. The insert shows such curves, taken at 
different temperatures. The different slopes before and after 
the transition temperature (dashed line at ^1.2 K) indicate a 
spin transition, due to a charge transfer onto the Gd3 system. 
The dotted line at very low temperatures (< 0.1 K) indicates 
a region where the J = 1/2 ground state of the Gd trimer, 
starts to become dominant (see text). Insert: Reversible, 
normalized magnetization curves measured at different tem- 
peratures, as indicated in the figure. 

with an abrupt change in slope at ~ 1.2 K. 

We fit the susceptibility data using a Heisenberg model 
for the Gd moments. The external magnetic field B cou- 
ples to the total moment J\ of Gd ion z, but the exchange 
interaction is mediated through the spin operators S\. 
Using S[ = (#i — 1) Ji, with g x the Lande g- factor of ion i, 
we can write the entire Hamiltonian as: 

H=~J2 ^4 ' $ + Ji J (9i ~ - 1) • Jh (1) 

i i<j 

where /iB is the Bohr magneton, Jjj is the exchange cou- 
pling between spins i and j = 1..3). We compute the 
susceptibility simply by enumerating all possible states 
using the quantum numbers J' = J\ + J 2 and J = J' + J3 . 
The thermal averaging of all eigenstates, with energies E^ 
and magnetizations Mk = dE^/dB, gives the susceptibil- 
ity as: 

X=^E M ^ Ek > ( 2 ) 
k 

where j3 — l/k&T, &b is Boltzmann's constant and Z = 
e~^ E]i is the partition function. 
In both the experiment and the fit, the susceptibility is 
normalized to the saturation magnetization at low tem- 
peratures (see the slope for T <1.2 K in Fig. 3). For the 



fit, the saturation magnetization from the three Gd ions 
is given by £)i=i...3 ^B9iy/jJ7i + 1). 

We find that the low-temperature susceptibility is well 
described by three Gd 3+ ions in 8 67/2 states {L = 0, S = 
7/2) interacting via equal exchange couplings = J « 
7.6 mK. For T > 1.2 K, it appears that one Gd changes to 
a 2+ state with a 7 Fq configuration {L = S = 3, J = 6). 
This shows that one electron is transferred from the Cso 
cage onto one of the Gd ions, leading to an overall spin 
transition for the Gd3 system. The model also agrees 
with a flattening of the 1/% vs. T curve at very low 
temperatures (below ~ 0.1 — 0.2 K) when the exchange 
couplings Jij becomes noticeable, and the J = 1/2 ground 
state of the system begins to dominate the susceptibility. 

In conclusion, we studied the magnetization dynamics 
of GdsN@C8o, a molecular system useful as a high con- 
trast agent for MRI, as a function of temperature. The 
magnetic studies are performed by using a sensitive on- 
chip SQUID technique. The initial susceptibility data 
are fitted with a quantum mechanical model, using all 
the states of the Gd trimer. At low temperatures, the 
data are explained by 3 Gd 3+ ions in state 8 5V/2 wea kly 
coupled antiferromagnetically. A spin transition of the 
Gd3 is observed at 1.2 K explained by a charge transfer 
of one electron from the Cso cage, onto one of the Gd 
ions. The later becomes a Gd 2+ ion, in a state 7 Fq. 
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